Hadrons and direct photon in pp and pA collisions at LHC and saturation effects 



O 



Amir H. Rczaeian^'^ and Andreas Schafer^ 

^ Institut fiir Theoretische Physik, Universitdt Regensburg, 93040 Regensburg, Germany 
Departamento de Fisica y Centra de Estudios Subatomicos, 
Universidad Tecnica Federico Santa Maria, Casilla 110-V, Valparaiso, Chile 

We investigate hadrons and direct photon production in pp and pA collisions at the energies of 
RHIC and LHC within the color-dipole approach employing various saturation models. We show 
that greatest sensitivity to saturation effects is reached at very forward rapidities for pp collisions at 
LHC (v^ = 14 TeV). The ratio of direct-photon to pion 7/71" production can be about 20 10 (at 
?7 = 7-^8 ). Therefore, direct photon production at forward rapidities should provide a rather clean 
probe. We calculate the rapidity dependence of the invariant cross-section and find some peculiar 
enhancement at forward rapidities which is more pronounced for direct photon production. We show 
that this peak is further enhanced by saturation effects. We provide predictions for the nuclear 
modification factor R^a for pions and direct photon production in pA collisions at LHC energy 
at midrapidity. We show within various saturation models that the pion Cronin enhancement at 
RHIC is replaced by a moderate suppression at LHC energy at midrapidity due to gluon shadowing 
effects. Cronin enhancement of direct photons can survive at LHC energy within models with a 
larger saturation scale. 
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I. INTRODUCTION 

The Large Hadron Collider (LHC) will allow to explore 
a new regime of QCD where parton saturation effects be- 
come important pJiSj]. At the same time, the physics of 
saturation might also be relevant for a detailed under- 
standing of the underlying events, i.e. the backgrounds 
for New Physics searches at LHC. 

It is believed that pp and pA collisions provide a test- 
ing ground to disentangle the initial- and final-state ef- 
fects in AA collisions and can be used as a baseline for 
understanding the physics of heavy-ions collisions. For 
example, to interpret jet-quenching, a precise and firm 
understanding of the Cronin, shadowing and saturation 
effects in pA collisions is indispensable. 

The Color Glass Condensate (saturation) approach to 
QCD at high energy [ll-lE| has been very successful to 
describe a variety of processes at Relativistic Heavy Ion 
Collider (RHIC) |6i](for a review see § and references 
therein). Nevertheless, the importance of saturation ef- 
fects is still disputable given that other approaches of- 
fered alternative descriptions, see for example Refs. 0,[1]. 
In order to test saturation physics and its relevance, it 
seems therefore mandatory to consider various reactions 
in different kinematic regions at LHC and future collider 
experiments. 

Here, we study hadron and direct photon production 
in pp and pA collisions within the light-cone color-dipole 
formulation and investigate the role of saturation and 
shadowing at LHC energies. The corresponding phe- 
nomenology is based on the universal qq dipole cross- 
section. The dipole cross-section incorporates the multi- 
ple gluon scattering and non-linear gluon recombination 
effects and can be in principle measured in deep-inelastic 
scattering (DIS), see section VI. In the parton model lan- 
guage, the dipole cross-section plays the role of leading 
twist parton distributions in an all twist environment. 

Direct photons (photons radiated in hadronic collisions 



not via hadronic decays) carry important information 
about the collision dynamics which is undisturbed by fi- 
nal state interactions. We compare hadron and direct 
photon production mechanisms at various energies and 
rapidities in pp collisions. We show that the ratio of 
photon/pion production at very forward rapidities grows 
and can become as big as one order of magnitude at 
the LHC energy ^/s = 14 TeV. Measurements of direct 
photons at forward rapidities should be rather clean, as 
the background from radiative hadronic decays is signif- 
icantly suppressed^. At the same time, we show that 
both hadrons and direct photons are sensitive to satu- 
ration effects at forward rapidities at ^/s = 14 TeV pp 
collisions. 

We also investigate the role of saturation and shadow- 
ing effects for hadron and direct photon production in 
pA collisions at LHC. Our approach gives a rather fair 
description of PHENIX data for the Cronin ratio RpA of 
pions. We show that the nuclear modification factor RpA 
for 7r° at LHC {\/s — 5.5 TeV) at midrapidity becomes 
less than 1 in all saturation color-dipole models due to 
gluon shadowing. The suppression obtained (for RpA) in 
our approach is less than the one predicted in the Color 
Glass Condensate (CGC) approach We will later 
highlight the difference between our results and other re- 
ported predictions. We will also show that the nuclear 
modification factor RpA for direct photons is also less 
than 1 within the CGC color-dipole model once shad- 
owing effects are included. In contrast, the Cronin en- 
hancement for photons can survive even after inclusion of 
shadowing effects within the Golec-Biernat and Wiisthoff 
color-dipole model which has a bigger saturation scale 
than the CGC model. 



^ Experimentally measurements at forward rapidities are a chal- 
lenge since production rates are lower due to kinematic limits. 
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The paper is organized as follows: In Sec. II and III 
we calculate gluon radiation from projectile gluons and 
quarks in the color-dipole approach. In Sec. IV we intro- 
duce the light-cone color-dipole factorization scheme for 
hadron production. In Sec III, IV we will also highlight 
the differences between our approach with others. In Sec 
V we calculate the direct-photon production in qN and 
pp{A) collisions. In Sec VI we introduce gluon satura- 
tion within various approaches and color-dipole models. 
In Sec VIII we discuss nuclear gluon shadowing, Cronin 
effect and nuclear modification factor for partons, pions 
and direct photon production. In Sec VII and VIII we 
present our numerical results for both hadron and direct- 
photon production in pp and pA collisions, respectively. 
As a conclusion, in Sec. IX we highlight our main results 
and predictions for LHC. 

II. GLUON RADIATION BY A PROJECTILE 
GLUON: gN{A) gigzX 

The underlying mechanisms of the multiple particle in- 
teractions is controlled by the coherence length Ic- In the 
incoherent case, the multiple interaction amplitude can 
be simplified as convolution of differential cross sections 
while in the coherent case, one should convolute scatter- 
ing amplitudes rather than differential cross-sections. 

The coherence length Ic can be estimated from the in- 
verse longitudinal momentum transfer, 

_ 2E, _ 2^,q(1 ~ a) 

where Ei is the initial parton energy and kx is the rel- 
ative transverse momentum of the final partons. In the 
above equation, M is the invariant mass of the two final 
partons, neglecting parton masses. The parameter a is 
the fractional light-cone momentum of one of the final 
partons. Gluon radiation is dominated by small values 
of a <C 1, therefore we have, 

Kj^ m nPt 

where Ef is the energy of the parton detected in the 
final state, px is the transverse momentum of the frag- 
mented hadron at midrapidity, and mAr is the nucleon 
mass. For pion production, the average momentum frac- 
tion (z) in the fragmentation functions is about 0.4 — 0.6 
in the range oi 2 < pT (GeV) < 8. For a coherence 
length which is shorter than the typical internucleon sep- 
aration Ic ^ Ra (where Ra denotes the nuclear radius), 
the projectile interacts incoherently. At the RHIC energy 
^/s = 200 GeV and intermediate pT we are almost in the 
transition region between the short- and long-coherence 
length regime. In more central collisions, at higher pT 
we are in the short-coherence length (SCL) limit and at 
LHC energies at moderate px we are again in the long- 
coherence length (LCL) limit. 




FIG. 1: Gluon radiation (top panel) and direct photon pro- 
duction (bottom panel) for a projectile quark interacting with 
the target. 



There is much experimental evidence for a large in- 
trinsic momentum of gluons, see Refs. [l^, [lH and ref- 
erence therein. Therefore, interaction with spectators is 
important since color screening is at work. At smaller 
and moderate px one should then include interaction 
with spectators, i.e. instead of "elastic" gluon scattering, 
gN — >■ gX, we need to consider bremsstrahlung subpro- 
cesses, gN — >• ggX, or qN — >• qgX. The lowest order 
for these processes includes the three graphs shown in 
Fig. [1] (interactions with the initial and two final par- 
tons). After summing over radiated gluons, the cross 
section of this reaction can be expressed in terms of the 
color-dipole amplitudes fio', 121 , and can be diagonalized 
for a nuclear target provided that the coherence length 
is sufficiently long. Note that since parton trajectories 
before and after gluon (or photon) radiation have differ- 
ent impact parameters, and the corresponding terms in 
the bremsstrahlung amplitude have different signs, one 
arrives at an expression, which is formally identical to 
the amplitude of an inelastic dipole-target interaction. 
This is only a formal procedure of calculation, while no 
real qq color-dipole is involved in the process of radiation 
contrary to DIS where a photon does split into a real qq 
pair. 

In the LCL regime, the transverse momentum spectra 
of gluon bremsstrahlung for a high energy gluon inter- 
acting with a nucleon N (or nucleus A) including the 
nonperturbative interactions of the radiated gluon reads 
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{kT,x) 



1 



X e 



(6,r2,a;)-A/-3^(^)(6,(fi-f2),x) 



(3) 



where a = P+{gi)/p+{g) ^ 1 denotes the hght-cone mo- 
mentum fractional of the radiated gluon. The partial 
amplitude N^g of a 3-gluons system colliding with a pro- 
ton at impact parameter b can be written in terms of the 
qq dipole amplitude [H, [3 , 



(4) 



where the factor 9/8 is the ratio of Casimir factors. Here 
the vectors f, ar and (1 — a)f denote the two gluon 
transverse separations r{gi) — f{g2), f{g) — r{g2) and 
r{g) — f{gi), respectively. 

Note that Eq. ([4]) can be simply understood by look- 
ing at several limiting cases: if r goes to zero, the 
transverse separation of final state gluons gi and 172 be- 
comes zero leading to M^g{b,r — 0,x) = which re- 
flects the fact that a point like gluon-gluon fluctuation 
cannot be resolved by interactions. In the two limit- 
ing cases of a — > 0,1, the three-gluon system will be 
reduced to the two-gluon system which can be then re- 
lated to the qq dipole cross-section via the Casimir factor: 
lima^o.i ■^i^{b, r, x) = jAf^{b, r, x) where 9/4 is the ra- 
tio of the octet and triplet color Casimir factor. 

The qq dipole amplitude in Eq. is related to the 
dipole-proton cross-section by integration over impact 
parameter. 



agg{r,x)=2 / d^bN^''Jb,f,x). 



(5) 



We still have to specify the light-cone distribution func- 
tion (^'gg) for the gg Fock component fluctuations of the 
incoming gluon, which includes nonperturbative interac- 
tions of these gluons. The light-cone wave function of 
the gluon-gluon (and quark-gluon) Fock component of a 
gluon (quark) was calculated in Ref. within a model 
describing the nonperturbative interaction of gluons via 
a phenomenological light-cone potential of an oscillatory 
form. This is given by. 



exp 



2r2 



K(e* •e)(e2* ■ r) 



+ (1 - a)(e* • e}{e* ■ r) - a(l - a){e* ■ e;){e- r) ^ 

where = 0.3 fm is the parameter characterizing the 
strength of the nonperturbative interaction which has 



been fltted to data on diffractive pp scattering [l^]. In 
Eq. IS]) the product of the wave functions is averaged 
over the initial gluon polarization, e, and summed over 
the final ones, ei_2. 

Based on pQCD one might expect that the gluon-gluon 
potential differs from the quark- ant iquark one simply by 
a Casimir factor 9/4. However, there exists plenty of ev- 
idence indicating that the interaction of gluons is much 
stronger due to non-trivial properties of the QCD vac- 
uum, see Ref. ^7] and references therein. It turns out 
that the exact shape of the light-cone gluon-gluon (quark- 
gluon) potential is not crucial [l3]- What is only impor- 
tant is the smallness of the mean quark-gluon separation 
ro which defines the effective strength of gluons interac- 
tion. The value of rn = 0.3 fm obtained from analysis 
of diffractive data [1^ agrees with both lattice calcula- 
tions [iHl and also with the phenomenological model of 
the instanton liquid [l6j . 

We consider here the asymptotic expression of the 
gluon radiation cross-section given in Eq. ([3]) for a —> 
which is reliable at very long coherence lengths. This is 
certainly valid at LHC energies. At RHIC energies, for 
hadrons produced at midrapidity with moderate pr, we 
are in the transition region between the regimes of long 
and short coherence lengths. Moreover, the color-dipole 
models we use in this paper, were fltted to DIS data at 
very small Bjorken-x xb < 0.01, which corresponds to 
£ 2 GeV at RHIC. Therefore, the prescription pre- 
sented here should be less reliable at high-p^ at RHIC 
energy. We will come back to this point in Sec. VIII. 

After some algebra one obtains. 



X 



9a, 



d^kr d^b 71-^ 



drUUb^r) 



X |J(l-e-4'^o/^)ji(fcT.) 
- Jo(fcTr)e^/(r)|, 

where the function /(r) is deflned as 



e 



(7) 



fir) = 



dA 



d9- 



(A2 -r2)Ar 



(A2+r2)2-4(Arcos(6l))2 



(8) 



In the case of a nuclear target the functional form of 
Eq. ([3|) still holds, but the dipole amplitude for a nu- 
cleon target should be replaced by the one for a nu- 
clear target A/"^. The partial elastic amplitude Af^ for 
a colorless three-gluon system colliding with a nucleus A 
can be written in terms of the partial amplitude Af^ of 
a three-gluon system colliding with a proton at impact 
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parameter b, 



cross-section is given by [10|, |lj 



^^qN{A)^qgX 



(9) 



where the 3-gluons amphtude A/"-^ is related to the qq 
dipole amphtude via Eq. ^ and 7/1(6) is the nuclear 
thickness function normalized to / (PbTA{b) = A. In a 
very similar fashion as for the nucleon target case, one 
can analytically carry out some of the integrals, 



(FkT(Pb V Jo \ kT\ 

X Ji(fcTr)e^"o +Jo(fcT'')e fv")) 



+ 



(2^ 

krp 



(10) 



with the notation. 



lG{b,r) = ^ / <fsN^^[s,r)TA{b + ^, 



aqg{r,x)TA(b), 



(11) 



where in the second line we used Eq. ([5]) and ignored pos- 
sible correlations between the color-dipole amplitude and 
nuclear thickness. Notice that the second line were iden- 
tically true if the nuclear profile would be a constant. 
In a more sophisticated approach in order to properly 
incorporate the correlation between the color-dipole am- 
plitude and the nuclear thickness, one should also have 
a model for the dipole amplitude which depends on the 
angle between the dipole transverse radius f and the im- 
pact parameter b. Unfortunately, with available HERA 
data, it is difficult to incorporate the color-dipole orienta- 
tion and most dipole models fitted to HERA data depend 
only on the absolute value of the transverse dipole size \r\ 
and impact parameter For a recent attempt to incor- 
porate the color dipole orientation, see Ref. UJ- It has 
been shown that the color-dipole orientation gives rise to 
azimuthal asymmetries |17l] , but is unimportant for total 
cross-sections. 

The remaining integrals in Eqs. (|7ll0p can be per- 
formed only numerically. 



III. GLUONS RADIATION BY A PROJECTILE 
QUARK: qN{A) qgX 

Gluon radiation of a projectile quark interacting with 
a nucleon(nucleus) qN{A) — )• qgX can be calculated in 
a similar way as outlined in the previous section. The 



{kT,x) 



1 



d^kT ^ " ' (2^) 
X e'^^-(^~^-^^^%;,(r-i,a)M',,(r2,a) 



(fb (fri (fr2 



^g<if^ ri,ri- af2,x) + Afgll^' {b, ^2 , ^2 - an , x) 



Ug^^\b, a{n ~ r2),x) - AAj^(-^)(&, in - r2), x) 



(12) 



where fi and r2 are the quark-gluon transverse sep- 
aration in the direct and complex conjugated ampli- 
tude respectively. For brevity, we define again a as the 
fractional LC momentum of the radiated gluon, a = 

p+i9)/p+{q) < 1- 

In Eqs. (|3I12I) we have already integrated over the 
transverse coordinates of the second parton. Note that 
the formulas in Eqs. pil2p are given in impact parame- 
ter representation and contain the sum of diagrams given 
in Fig. [TJ The derivation of these equations can be found 
in Refs. [l3, [13 • The collinear divergences which are 
the source of scale dependence of the parton distribution 
functions and fragmentation functions in the factoriza- 
tion Eq. (PO)) are already subtracted in these equations. 

The interaction amplitude of a colorless gqq and gg 
system with a nucleon target can be written in terms of 
qq dipole amplitudes [isj . 

A/;^,(6,rl,r2,a;) = ^{a/-^(6, fi, 2;) + A/-^(&, r^, x)} 



AAj^(6,f,a;) 



■A/-i^(6,ri-r2,x), 
Mb,f,x). 



(13) 
(14) 



Again likewise Eq. (|4]), the above equations immediately 
satisfy several simple limiting cases. When the qq trans- 
verse separation goes to zero i.e. ri ~ ^2, the qq pair is 
indistinguishable from a gluon, and Eq. (1131) correctly re- 
duces to N^qqib, rx,rx,x) — jj\f^{b, fi,x). Moreover, in 
the limit of vanishing fi (or r2), the qg (or qg) system is 
indistinguishable from a quark (antiquark) and Eq. (|13p 
becomes Af^^^ {b,ri,0,x) ^ ' " 

In the derivation of Eqs. (131 [T^ . one can rearrange 
the final result in terms of Eqs. (|4I13I14I) . However, 
this is more than just some change of notation since the 
combination of the right-hand sides of Eqs. (|4ll3ll4p are 
already well-known as forward scattering amplitudes of 
ggg^ gqq and gg system interacting with a proton tar- 
get. Eqs. (j4ll3ll4p are exact and are not based on any 
approximation. A formal derivation of these equations 
is similar to the derivation of the qq-proton dipole cross- 
section, namely one replaces the qq system by a qqg (or 
999 I 99) 0116. The exchanged gluons can now couple to 
different partons in the qqg system (or ggg , gg) which 
generates different phase factors. The precise calcula- 
tion of the color traces for the different couplings of the 
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exchanged gluons to the quark, antiquark and gluon ( 
or 3-gluons, gluon-gluon) leads to the exact expression 
given in Eqs. (|4ll3ll4p ^ [Tj-ll l Notic e that in the 
CGC approach the relations Eqs. (|4ll3ll4p holds only if 
one assumes that the weight function for averaging over 
the target color charges are Gaussian [l8l |. 

The forward scattering amplitude of qq, gqq and gg 
interacting with a nucleus target at impact parameter 6, 
can be again written, in eikonal form, in terms of the 
dipole elastic amplitude J\f^ of a qq dipole colliding with 

a proton at impact parameter &, 

MM,r,x) - i_e-/^'^^^i^(^»)^-(^+^"^, (15) 



"qq 

J^^qqib,r,x) 



J\f^g{b,r,x) = i_e-3/'i^^^A^f,(^»)r.(ft+«-). (17) 

The light-cone distribution of quark-gluon fluctuations 
^'gg in Eq. (IT^ is given in Ref. 10]. In the hmit a <C 1 
which is of practical interest at high energy, the quark- 
gluon distribution function including non-perturbative 
effects has the form, 



99 



(f,a) = --J^— exp(-rV2r2). (18) 



where the parameter tq = 0.3 fm denotes the mean 
quark-gluon separation and is the result of a fit to soft 
diffraction pp — > pX. 

One can show that for a <C 1 the cross-section of 
gluon bremsstrahlung from projectile quarks is 6 times 
smaller than the corresponding cross-section for a pro- 
jectile gluon given by Eqs. p|6ll2ll8p due to the color 
factor: 



^qN{A}-^qgX 



^gN{A)^gig2X 

6 ■ 



(19) 



Note that similar results as Eqs. (131121) w as also ob- 
tained by Jalian-Marian and Kovchegov [l9| in a color 
glass condensate picture where the color dipole ampli- 
tudes in Eqs. (|3ll2p are replaced by a product of two 
Wilson lines evaluated in the field of the color glass con- 
densate. See also Ref. [2^ for an earlier attempt along 
this line. Loosely speaking, these two formulations are 
equivalent in the quasi-classical (Glauber) approxima- 
tion. However, in order to include small-x evolution, it is 
not sufficient to only put Wilson lines in the evolved CGC 
fields. This only leads to logs of energy in the rapidity in- 
terval between the produced gluon and the nucleus. One 
should also include the evolution in the rapidity interval 
between the projectile and the produced gluon [2l| , thus 
describing gluon emission. It was shown by Kovchegov 



and Tuchin that such an evolution is the linear BFKL 
equation due to some very interesting cancellations of all 
nonlinearities. 

In our approach, the effects of gluon emissions between 
the quark (gluon) and the produced gluon (and its evo- 
lution) are effectively included in the master Eqs. (I3ll2p 
via the non-perturbative quark-gluon (gluon-gluon) light- 
cone distribution functions Eqs. (|6ll8p which is obtained 
from a fit to soft pp diffraction data. The diffractive 
excitation of the incident hadrons to the states of large 
mass is a more sensitive probe of gluon-gluon fluctuations 
than the total cross section [l^]- While the gluon emis- 
sions between the projectile and target including their 
non-linear recombination effects are effectively incorpo- 
rated in terms of color-dipole forward amplitudes ob- 
tained from a fit to DIS data. By means of Eqs. pil2p 
one can also describe the long-standing problem of the 
small size of the triple-pomeron coupling [lol . [ll| . 

A word of caution is in order here. Notice that al- 
though the non-perturbative gg and qg light-cone distri- 
bution functions include some saturation effects of the 
projectile proton UII- Nevertheless, the gluon pro- 
duction cross-section given by Eqs. (|3I12I) is intrinsically 
asymmetric, namely it treats the "projectile" proton ap- 
proximately in a collinear factorization framework while 
treating the "target" proton (or nucleus) in a saturation 
framework. Strictly speaking this may be justified only 
in the case when saturation effects are present in the tar- 
get wave function, but are absent in the projectile wave 
function, such as in pA collisions or in forward particle 
productions. Although it appears that such a simple ap- 
proximation is sufficient to describe the existing exper- 
imental data for hadron and direct photon production 
at small x at midrapidity in pp collisions, see Figs. [Sj |6l 
Nevertheless, our formulation at midrapidity in pp col- 
lisions is not well justified and therefore our results at 
midrapidity in pp collisions may not be valid. 



IV. HADRONS PRODUCTION IN 
HIGH-ENERGY pp AND pA COLLISIONS 

The cross section of hadron production in pp (or pA) 
collisions at impact parameter b is given by a convolu- 
tion of the distribution function of the projectile gluon 
or quark inside the proton with the gluon radiation cross- 
section coming from gN or qN (gA or qA) collisions and 
also with the fragmentation functions. For simplicity, we 
assume here that the projectile gluon/quark has the same 
impact parameter relative to the target as the beam pro- 
ton. This is certainly a rather poor approximation which 
we will try to improve upon in future. 
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r-1 



-■ - z a^kgTao z z z^ 

(20) 



fq/p{xq, Q^) and fg^p{xg, Q^) are the parton distribution 
functions (PDF) of the colHding protons, which depend 
on the hard scale Q and the hght-cone momentum frac- 
tions Xq and Xg for quarks and gluons, respectively. The 
function Dfi/q g{z,Q'^) is the fragmentation function of 
parton q, g to the final hadron h with a momentum frac- 
tion z. In the above equation, the variables ^ and ^ are 
momentum fractions of a parton in the beam and target. 
The variables Xi^2 are defined by. 



Xi 



X2 = ^e-\ (21) 



where pT and rj are the transverse momentum and rapid- 
ity of the produced hadron. 

In Eq. (PH)) the cross-sections of gluon radiation in 
gp{A) — 7^ ggX and qp{A) — qgX are given by Eqs. pil2p . 
We assume that the projectile parton acquires high 
transverse momentum fcy as a result of coherent mul- 
tiple rescattering, while the radiated gluons that gener- 
ate this momentum are summed to build up the color 
dipole cross-section. Then, explicit inclusion of gluon 
bremsstrahlung balances the large kx- 

Notice that in the dipole approach in contrast to the 
parton model, one should rely on the parton distribution 
functions taken at a soft scale since the evolution to the 
hard scale is performed via gluon radiation, which is en- 
coded in the phenomenological dipole cross-section fitted 
to DIS data for the proton structure function. However, 
the dipole cross-section misses the Q^-evolution of the 
xi -distribution, which is especially important at forward 
rapidities, since the parton distributions fall off at a;i — > 1 
much steeper at high . In order to account for this ef- 
fect and provide the correct xi-distribution, we take the 
integrated parton distribution in Eq. (j20p at the hard 
scale Q = kr 22, 23]. 

Notice that at high energies and midrapidity the par- 
ton fractional momenta in the beam and target are small, 
Xi ~ a;2 <C 1, so hadron production is dominated by frag- 
mentation of radiated gluons gp{A) gig2X. However, 
at very forward rapidities the quark contributions are 
important and the subprocess qp{A) — ^ qgX becomes 
relevant. Therefore, different subprocesses dominate in 
different kinematic regimes and their overlap is small. 



V. PHOTON RADIATION IN HIGH-ENERGY 
pp AND pA COLLISIONS 

Production of direct photons in the target rest frame 
should be treated as electromagnetic bremsstrahlung by 
a quark interacting with the target. In the light-cone 
dipole approach the transverse momentum distribution 
of photon bremsstrahlung by a quark propagating and 
interacting with a target nucleon (or nucleus A) at impact 
parameter b, as calculated from the diagrams in Fig. [T] 
(we show only the single gluon exchange diagrams), can 
be written in the factorized form [13, u% u31 



daiqNiA)^jX) 
d{lna)d^PTd'^b (2 



7r)2 51 / 



X e 



ri)0^,(a, r2) [^fqf^\b, an,x) 

+ Afqf^\b,af2,x)-N-q''q^^\b,ain - r2),x)] , (22) 

where fi and r2 are the quark-photon transverse sepa- 
rations in the direct and complex conjugated amplitudes 
respectively; a — p^ /p^ denotes the fractional light-cone 
(LC) momentum of the radiated photon. Correspond- 
ingly, the transverse displacements of the recoil quarks 
in the two amplitudes are ari and ar2 respectively. In 
Eq. (P^. (fi^qia, f) is the light-cone (LC) distribution am- 
plitude of the projectile quark 75 fluctuation. Averaging 
over the initial quark polarizations and summing over all 
final polarization states of the quark and photon, we get 

^(I)*g{a,ri)(j)jq{a,r2) = a^Ko{amqri) 

inj I 

X Kf)(am.nr2) + \l + {I — a)^]— — -Ki(amqri) 

7-1 f 2 



X Ki{amqr2) 



(23) 



where Ko,i{x) denotes modified Bessel functions of the 
second kind and niq is an effective quark mass, which 
can be regarded as a cutoff regularization. Following 



Refs. 



we take m„ = 0.2 GeV. The forward scat- 



tering amplitude Af^ can be again written, in the eikonal 
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form, in terms of the dipole elastic amplitude of a 

qq dipole colliding with a proton at impact parameter b 
as defined in Eq. ([T5|) . 

In order to obtain the hadron cross-section from the el- 
ementary partonic cross section Eq. (I22[) . one should sum 
the contributions from quarks and antiquarks (since only 
quarks and antiquarks can radiate photons) weighted 
with the corresponding parton distribution functions. 
The PDFs of the projectile enter in a combination which 
can be written in terms of proton structure function 
F2{x,Q^). Notice that the contribution of gluon split- 
ting to quark- antiquark pairs (and higher Fock compo- 
nents) is already contained in the sea quark distributions 
of the proton. Therefore, the direct-photon production 
cross-section in pp and pA collisions is given by [23 - l23 | , 

d(j{pp{A) jX) _ xi da 
dxpd'^PT d'^b xi + X2 J XI oP- 



^zKg/(^) + g7(-)} 



xi d(7{qp{A) ^ ^X) 
a d{lna)d'^pTd% 



{PT,X2), 



1 



Xi + X2 



da{qp{A) -> -yX) 
d{lna)d^PTd^b 



{■PT,X2), 

(24) 



where the variable xi and X2 are defined in Eq. (j2ip 
and xp = Xi — X2 is the Feynman variable. We have 
recently shown that in this framework one can obtain 
a good description of the cross-section for prompt pho- 
ton production in pro ton-proton collisions at RHIC and 
Tevatron energies |22l . [23j . and Drell-Yan dilepton pair 
production [22l [25j. Here, we employ this formulation to 
give predictions for the ratio of photon/pion production 
cross-sections at various rapidities for LHC. We will also 
provide prediction for the nuclear modification factor in 
pA collisions at LHC. 

Notice that in the color-dipole factorization 
Eqs. (I20I24P neither K-factors (next-to-leading-order 
corrections), nor higher twist corrections should be 
added. The phenomenological dipole cross-section fitted 
to DIS data should already incorporate all perturbative 
and non-perturbative radiation processes. The only 
contribution which is still missing in Eq. (|24p is the 
effect of the primordial momentum of the projectile 
parton. However, it has been shown that in the color- 
dipole approach, the primordial momentum should have 
a purely non-perturbative origin, and is considerably 
smaller than in the parton model [12, [2^. This effect 
should be of little importance for the kinematic regions 
of interest of this paper. 

A word of caution is in order here. The type of fac- 
torization scheme outlined above Eqs. (|20l24p has not 
been yet rigorously proven at any order of pQCD in the 
kinematic region of our interest and is most probable 
not exact. Nevertheless, there is growing evidence in the 
literature that it gives a good approximation for the pro- 
cesses discussed here [3, US [13, l23 - l28| . 



VI. GLUON SATURATION AND COLOR 
DIPOLE MODELS 



At high energies /small Bjorken-x, QCD predicts that 
gluons in a hadron wavefunction form a new state, the 
so-called Color Glass Condensate (CGC) il-5]. The cor- 
nerstone of the CGC is the existence of a hard saturation 
scale Qs at which nonlinear gluon recombination effects 
become important and start to balance gluon radiation. 

The concept of saturation and the taming of the power- 
like rise of the gluon distribution at small x was first ad- 
dressed by Gribov, Levin and Ryskin in the double loga- 
rithmic approximation [ij . A first hint toward saturation 
effects at HERA came from the phenomenologically suc- 
cess of the Golec-Biernat and Wiisthoff (GBW) model 
[29} . This model incorporates the basic saturation effects 
into the color-dipole cross-section on a proton target. In 
the CGC framework the dipole-proton forward scattering 
amplitude can be in principle found by solving the per- 
turbative nonlinear small-x Balitsky-Kovchegov (BK) [^ 
or Jalilian-Marian-Iancu-McLerran-Weigert-Leonidov- 
Kovner (JIMWLK) [3] quantum evolution equations. 
The BK and JIMWLK evolution equations unitarize the 
linear Balitsky-Fadin-Kuraev-Lipatov (BFKL) [131 evo- 
lution equation at small-x in the large-A'c limit (BK) 
and beyond (JIMWLK). It has been shown that next-to- 
leading-order (NLO) corrections to the BFKL equation 
(and therefore to BK and JIMWLK kernels) are large 
and negative [3l|. There was no reason to believe that 
still higher order corrections are unimportant, until quite 
recently, when it was found that the consistent incorpo- 
ration of the running coupling Ug into the BFKL, BK and 
JIMWLK equations p32H34| leads to phenomenologically 
rather successful descriptions. Still the actual calcula- 
tion of higher-order corrections to these non-linear evo- 
lution equations remains as a challenge. Thus, we resort 
to a QCD-like model which incorporates the basic fea- 
tures of gluon saturation into the dipole-proton forward 
scattering amplitude, and provides predictions which will 
allow to test the validity of our treatment. There are sev- 
eral parametrizations proposed in the literature which all 
give a good description of HERA data but predict differ- 
ent saturation scales, see Fig. [2] In this section we re- 
view some of these models and later we will employ them 
for hadron and photon production in various kinematic 
regimes and investigate the uncertainties of the various 
models and discuss the differences between them. 



A. GBW model 

The dipole-proton cross-section aqq{r,x) is usually 
written as an integral of the imaginary part of the for- 
ward scattering amplitude Af^{'r,b, s) over the impact 
parameter b as defined via Eq. ([5]). One may neglect 
the 6-dependence in Af^ making the integral in Eq. ([S]) 
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trivial, giving the proton's transverse area factor: 

aqg{r,x) EE aoAf^^{r,x). (25) 

A popular parametrization for the qq dipole cross-section 
on a nucleon target is due to Golec-Biernat and Wiisthoff 
(GBW) and is able to describe DIS data with a sim- 
ple form for the color dipole amplitude, 



GBW 



(26) 



where the x-dependence of the saturation scale is given 

by 



(27) 



The main feature of the model is that for decreasing x, 
the dipole amplitude saturates at smaller dipole sizes. 
Note that there is no unique definition for the saturation 
scale in literature. Following Refs. [1^, Issl - fSTj we define 
the saturation scale Ql = as a energy scale at which 
the qq dipole scattering amplitude Af becomes sizable. 



Mgg{r, = V2/Q„ x) = l- « 0.4. 



(28) 



For the GBW model, this definition coincides with the 
saturation scale Qs defined in Eq. (P7)) . The value of the 
intercept A ~ 0.25 — 0.30 is consistent with perturbative 
predictions based on small-x evolution [s^, 113, [38l - l40| . 
The parameters ctq = 23.9 mb, xq = 1.11 x 10~^, and 
A = 0.287 were determined from a fit to F2 for x < 
0.01 and G [0.25,45] in the presence of charm quarks 
with mass rric = 1.4 GeV [s^. Note that the saturation 
scale in the GBW model reduces with the inclusion of 
the charm quark [36j . 



B. CGC, AAMS-BK and b-CGC models 
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The linear DGLAP evolution equation which only in- 
cludes gluon radiation may not be appropriate for the 
saturation regime where nonlinear recombination subpro- 
cess are important. lancu, Itakura and Munier proposed 
an alternative color glass condensate (CGC) model (4l| . 
based on the BK equation . In this model the qq dipole 
amplitude for a nucleon target is parametrized as. 



CGC 



{r,x) 



2(7^ + 7^^ In ;4j) 



rQs < 2 

-Aln-iBrQ^) . ^Q^ > 3 ' 

(29) 

where the saturation scale is again parametrized as 
Eq. (HH), Y = ln(l/a;), and k = x"{ls)/x'{ls) where 
X is the LO BFKL characteristic function. The coeffi- 
cients A and B in the second line of are determined 
uniquely from the condition that the color dipole cross- 
section and its derivative with respect to rQs are contin- 
uous at rQs = 2: 



A 



B 



(l-TVo) 



(30) 



FIG. 2: Top panel: Saturation scale defined via Eq. (|28l) as a 
function of 1/x for various color-dipole models. Lower panel: 
The total dipole- proton cross section Gqq[r,x) at fixed x = 
lO"'^ in the various color-dipole models introduced in Sec. 
VI. 



The parameters 7^ = 0.63 and k = 9.9 are fixed at the 
LO BFKL values. The others parameters Ao = 0.7, ctq = 
35.7 mb, xo = 2.7 x 10""^ and A = 0.177 were fitted to F2 
for X < 0.01 and Q^ < 45 GeV^ and including a charm 
quark with nic = 1.4 GeV. Notice that for small rQs < 2, 
the effective anomalous dimension 1 — 7^ in the exponent 
in the upper line of Eq. (I^Hl) rises from the LO BFKL 
value towards the DGLAP value. 

Recently, Albacete, Armesto, Milhano and Salgado cal- 
culated numerically the dipole-proton scattering ampli- 
tude from the BK equation including running coupling 
corrections (AAMS-BKl, 2 model) |3j|. Note that the in- 
corporating of the running coupling is essential in this ap- 
proach, though its implementation is model dependent. 
The free parameters in their fit to HERA data are related 
to the initial condition for the evolution at x.;„ = 10~^. 
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They used two families of initial conditions, the GBW 
form (AAMS-BKl model) 



1 — exp 



(31) 



and the McLerran-Venugopalan form (AAMS-BK2 
model) : 



A/'„ ^(r,Xi„) = 1-exp 



In 



1 



■A 



QCD 



(32) 

where Q^g the initial saturation scale. In their global 
analysis of HERA data there are four free parameters 
which are fitted to i^2-data for x < 0.01 and Q^/GcV^ G 
[0.045,800]: the initial saturation scale Qso, the overall 
normalization (Tq: the infrared parameter C introduced 
in the running coupling and the anomalous dimension 
7. The values of parameters can be found in table 1 of 
Ref. dl. 

The gluon density is larger in the center of a proton 
b — than at periphery 6^2 — 3 GeV~^ probed in 
the total 7*p cross-section. Therefore, impact-parameter 
dependence of the dipole-proton forward scattering am- 
plitude seems to be essential. There has been several 
attempts to model the impact-parameter dependence in 
dipole-proton forward scattering amplitudes. We con- 
sider here the model proposed by Watt and Kowalski (b- 
CGC) [42]. In this model, the dipole-proton forward scat- 
tering amplitude has the same form as the CGC model 
Eq. (P^ . but the saturation scale Qs now depends on 
impact parameter. 



Qs 



Qs{x,b) = ( 



X / 



exp 



2B, 



CGC 



(33) 



The parameter i?cG = 7.5GeV is fitted to the t- 



dependence of exclusive J/^ photoproduction. It has 
been shown that if one allows the parameter 7s to vary 
together with the other parameters (in contrast to the 
CGC fitting procedure where 7s is fixed to its LO BFKL 
value) , this results in a significantly better description of 
data for F2 with the value of 7s — 0.46, which is remark- 
ably close to the value of 7s = 0.44 recently obtained 
from the BK equation Other parameters obtained 

from the fit are: Wq = 0.558, a;o = 1-84 x 10"^ and 
A = 0.119 4.42]. 

Notice that calculation of the pT-distribution of pro- 
duced hadrons/photons in pp collisions needs only knowl- 
edge of the total dipole cross-section and is independent 
of the impact-parameter dependence of the forward scat- 
tering dipole-proton amplitude. Nevertheless, the inte- 
grated dipole cross-section of the b-CGC model is differ- 
ent from other dipole models 



C. KLR-AdS/CFT model 

The above mentioned dipoles models are motivated by 
pQCD and their validity at very small where one has 
to consider small-x evolution in the large coupling limit 
is questionable. Performing calculations in the strong 
coupling limit of QCD is very difficult. One may re- 
sort to other QGD-like theories, such as A/" = 4 Super- 
Yang-Mills where one can perform calculations in the 
non-perturbative limit of large 't Hooft coupling by em- 
ploying the Anti-de Sitter space/conformal field theory 
(AdS/CFT) correspondence 14411. On this line, recently, 
Kovchegov, Lu and Rezaeian |45| proposed a new color 
dipole parametrization inspired by the AdS/CFT ap- 
proach (KLR-AdS/CFT) which reasonably well describes 
the HERA data for inclusive structure functions at small- 
X and Q^. In this model, the dipole-proton scattering 
amplitude is given by. 



Af,fHr,x) = 1-exp 



(- + - 



A{)Xr 



'l-x 



(34) 



with notations 



A = 



3mA 



A 



1 



1 



m< 2^ 



m > 27 



2m V 4m2 27nfi 
Ml{l~xY 



1/3 



where ^0 = V Xym GeV. The parameters of the model 



for quark mass 



140 MeV and 't Hooft cou- 



pling Xym ~ 10 obtained from the fit to the HERA 
data (in the range of a; e [6.2 x 10^^, 6 x 10^^] and 
gVCeV^ e [0.045,2.5]) are: Mq = 8.16 x lO^^ and 
(To = 26.08 mb (see Eq. (1^51) ). We will also consider an- 
other fit to the same data but with 't Hooft coupling 
Xym — 20 which also gives a good fit: Mq — 6.54 x 10^^ 
and (To = 22.47 mb [45l |. 



cos(^) 



27m 



(35) 



Similarly, the saturation scale in the KLR-AdS/CFT 
dipole model ([M)) can be obtained from the definition 
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given in Eq. ((281 
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-lAdS 



(x) 



2^0 2: 



A^5(l-x)7r \p: 



2Mo 

Prn 



l-x 



(36) 

In this model the saturation scale varies in the range of 
1 3 GeV becoming independent of energy/Bjorken-x at 
very small x (see Fig. [2]). This leads to the prediction of 
x-independence of the F2 structure function at very small 
X and in a region where there is no experimental data 
yet. 

Note that the KLR-AdS/CFT dipole scattering_^am- 
plitude exhibits the property of geometric scaling [46j : it 
is a function of rQ^'^'^[x) only. Moreover, the anoma- 
lous dimension in this model is 7^ = 0.5 which is rather 
close to the value of 0.44 obtained from the numerical 
solution of the BK equation ^43|]. Thus in many ways 
the predictions of the KLR-AdS/CFT model are simi- 
lar to the predictions of the CGC model. Therefore, the 
non-perturbative KLR-AdS/CFT model which is valid at 
low < 2.5 GeV^ could be viewed as complementary 
to the perturbative description of data based on satura- 
tion/Color Glass Condensate physics. The main differ- 
ence is the x-dependence of the saturation scale Q^'^^{x), 
which leads to x-scaling at small x and . 



D. Semi-Sat Model 
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FIG. 3: The hadrons spectra at the RHIC energy at midra- 
pidity. The dashed lines are calculated with the GBW model. 
The solid lines are the pQCD calculation results taken from 
Ref. [HH]. Note that the color dipole approach is valid at 
very small X2 corresponding to pr < 2 GeV at RHIC energy 
and midrapidity (shown by a line). The experimental data are 
from [53, ihM ■ 



In order to demonstrate the importance of saturation, 
we will also use a semi-saturation model (Semi-Sat) fitted 
to F2 with X < 0.01 and £ [0.25,45] GeV^: 

_^Scmi-Sat(^-.^ fe, x) = 2AAo (^^) ^^^^ , (37) 

where Qs is defined in Eq. ([33l) . The parameter jeff 
is defined for rQs < 2 a.s -fef f = Is + In^, and 
for rQs > 2 as Jeff — 7s- The other parameters are 
given by 7^ = 0.43, Wq = 0.568, xq = 1.34 x lO^^ and 
A = 0.109 42]. Surprisingly, the fit obtained with such 
an oversimplified model is as good as for the other models 
with xVd-o.f. = 0.92. 

Comparing Eq. ([29]) and Eq. (|37|) one can see that 
they treat the region rQg > 1 differently. The CGC 
model describes this region based on solutions to the BK 
equation [47l - l49| for rQs>2 (with a phenomenological 
matching at r Qs = 2) which are also applied (somewhat 
inconsistent) in this model for r close to 1/Qs. 

In Fig. [21 we show the saturation scale (top panel) and 
qq dipole-proton cross-section (lower panel) within vari- 
ous color-dipole models fitted to the HERA data. Note 
that we used for all curves in Fig. [5] the same definition 
for the saturation scale given in Eq. ([25)) . It is obvious 
that the discrepancies among different models fitted to 
the same data are quite significant. Therefore, it seems 
that HERA data alone is not sufficient for a satisfactory 
understanding of saturation physics. One of the aims of 



this paper is to investigate if hadrons and photon produc- 
tion at LHC can improve our understanding of saturation 
effects. 



VII. NUMERICAL RESULTS FOR pp 
COLLISIONS 

In order to analytically reduce the four-dimensional 
integrals in the partonic cross-sections Eqs. pil2p to 
one-dimensional integrals Eqs. ()7ll0p . we assumed that 
the strong coupling as is a constant. In principle, the 
strong coupling as entering in the gg and gq light-cone 
distribution functions of the incoming parton defined in 
Eqs. (I6|18|) is a function of the transverse dipole size. To 
improve our description, we replace as by asikx), where 
kx is the transverse momentum of the parton. More pre- 
cisely, in Eqs. ()3I12I) we replace as(ri)as(r2) — a'^{kT) 
where ri and r2 are the gluon-gluon (or quark-gluon) 
transverse separation in the direct and complex conju- 
gated amplitudes respectively and are related by a dou- 
ble Fourier transformation to the transverse momentum 
of the radiated gluon fc^, see Eqs. (j3|12|) . 

We employ recent NLO parton distribution functions 
(PDFs) developed for LHC apphcation (MSTW2008) 
[50|. For the fragmentation functions (EEs) we use the 
result of a recent NLO AKK08 analysis [51[. For the 
running strong coupling Og, we employ the same scheme 
as used for the MSTW2008 PDFs, namely we solve the 
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FIG. 4: Pions and protons spectra in pp collisions at the 
LHC energies ^/s — 5.5 and 14 TeV at midrapidity. Theory 
curves are calculated with the GBW model. Note that our 
model for particle productions is not reliable at midrapidity, 
for explanation see the end of Sec. III. 



renormalization group equation in the MSbar scheme at 
NLO level [11] • We stress that all phenomenological pa- 
rameters in our model are already fixed by other reac- 
tions and in this sense our results can be considered as 
parameter-free predictions. 

In Fig. [3J we show dipole model results obtained 
from the light-cone factorization in Eq. for pion 

(Tr^jTr"*" + TT~) and proton {p + p) spectra at RHIC en- 
ergy y/s = 200 GeV and midrapidity. The experimen- 
tal data are from PHENIX [H and STAR [5|. For a 
compassion, we also show the results coming from an 
improved pQCD calculation performed in Ref. [s^. No- 
tice that in the parton model results shown in Fig. |3] a 
fixed A'- factor K = 1.5 was introduced in order to sim- 
ulate higher order perturbative corrections while in the 
color dipole approach we do not introduce a if-factor 
since the dipole-proton cross-section fitted to HERA in- 
corporates all higher order radiations. Note that all the 
above-mentioned parametrizations for the color dipole 
cross-section have been fitted to DIS data a.t x < 0.01. 
This corresponds to pt < 2 GeV for RHIC energy at 
midrapidity (see Eq. ©), so the PHENIX and STAR 
data plotted in Fig. |3]are not suited for a model test. It 
is seen from Fig. |3] that deviation of color dipole results 
from the experimental data starts at about pt = 2 — 4 
GeV. At LHC energies = 5.5 and 14 TeV for a 
large range of pr (even at 77 = 0) we have X2 ^ 0.01, 
therefore we expect the color dipole prescription to be 
valid. In Fig. 21 we show the predictions of the GBW 
model for pion spectra in pp collisions for LHC energies 
^/s = 5.5, 14 TeV at midrapidity 77 = 0. The predictions 
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FIG. 5: Pion spectra obtained from various dipole models at 
forward rapidities and LHC energies for pp collisions. 

for pion invariant cross-sections at various rapidities in 
pp collisions for LHC are given in Fig. [5l One can see 
from Fig. [5] that various dipole models presented in the 
previous section with explicit saturation give rather sim- 
ilar results (we will scrutinize this below). Note that the 
KLR-AdS/CFT model described in Sec. VLC was fit- 
ted to the HERA data with x G [6.2 x 10""^, 6 x 10"^] 
and QVGeV^ G [0.045, 2.5]. Therefore it is only vahd at 
very forward rapidities and low pT- As it is seen in the 
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FIG. 6: Direct photon spectra obtained from the GBW dipole 
model at the RHIC and CDF energies for pp collisions. We 
also show the NLO pQCD curve from the authors of reference 
[111 (given in table 3 of Ref. ^) which used the CTEQ5M 
parton distribution functions with all scales set topT- Exper- 
imental data are from the PHENIX exp eriment [Sa ] at 77 = 0, 
and from the CDF experiment [13, at I77I < 0.9. The er- 
ror bars are the linear sum of the statistical and systematic 
uncertainties. 



GBW 
CGC 

AAMS-BKl 

KLR-AdS/CFT, X,.,. =10 



VS"= 14 TeV, pp^ic" X 



1 GeV 




FIG. 7: Invariant cross-section for pion (top) and direct pho- 
ton (lower panel) production in pp collisions at LHC as a func- 
tion of rapidity 77 calculated with various color dipole models 
for various fixed pr- Our results for pp collisions is less reliable 
at midrapidity (shown with brown color). 



upper panel of Fig. [SJ the two color dipole solutions of 
the BK equation for the GBW and MV initial conditions 
(AAMS-BKl, 2) give very similar results and further on 
we will only consider one of them. 

In Fig. [51 we show direct photon spectra obtained in 
our color-dipole approach Eq. (l24l) . at the RHIC ^] 
= 200 GeV) and GDF (Vi = 1-8 TeV) energy 
[stI [59I I . Again, we should warn that our results at high 
Pt for lower energies like RHIC and CDF are less reli- 
able since X2 > 0.01 which is beyond the limit of applica- 
bility of the color-dipole light-cone factorization scheme. 
Nevertheless, the agreement of our results with available 
data for both hadron and photon production at RHIC 
and CDF energies is rather satisfactory for x < 0.01. As 
a comparison, in Fig. [61 we also show the NLO pQCD 
curve for CDF energy [56|. The predictions for direct 
photon spectra at LHC energies in pp collisions within 



various color-dipole models can be found in Ref. |23l |. 

In Fig. [3 the differential cross-section of pion 7r° (top 
panel) and direct photon 7 (lower panel) production at 
LHC are plotted versus rapidity at fixed transverse mo- 
menta Pt = 1 and 2 GeV within various color-dipole 
models. It is seen that the discrepancies among vari- 
ous saturation color dipole model results can be about 
a factor of 2 — 3 at moderate rapidities. At the kine- 
matic limit, i.e. at very forward rapidities and higher 
Pt where the differential cross-section approaches zero, 
kinematic constraints limit the parton phase space and 
saturation effects become less important. This is seen 
in Fig. [7] where as we approach very forward rapidities 
at the kinematic limit, the discrepancies among various 
saturation models shrink, and the invariant cross-section 
identically approaches zero. Notice that for hadron pro- 
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FIG. 8: Invariant cross-section for pion production in pp 
collisions at LHC as a function of rapidity rj calculated with 
the GBW color dipole model for various fixed pr- 



FIG. 9: Invariant cross-section for direct photon production 
in pp collisions at LHC as a function of rapidity r) calculated 
with the GBW color dipole model for various fixed pt- 



duction in the master Eq. (I^Ul) . the hght-cone momen- 
tum fraction x = ^ (where < z < 1 is the fragmenta- 
tion fraction) enters the gluon radiation cross-section and 
therefore the color dipole cross-section, while in the case 
of direct photon production Eq. (j24p . we have x = X2- 
Therefore, the applicability of the KLR-AdS/CFT model 
which is valid for x < 6 x 10""'^ (and < 2.5 GeV'^), 
can be extended for direct photon production to lower 
rapidities compared to the case of hadrons. It is seen 
from Fig. [7] that for both hadron and photon produc- 
tion, away from the kinematic limit, at not very large 77 
and Pt, a color-dipole model with larger saturation scale 
leads to a stronger peak at forward rapidity (having in 
mind that the saturation scale is a dynamical function of 
X, see Fig. [2]). 

In Fig.[7J it is seen a peculiar enhancement of the pho- 
ton production rate at forward rapidities. This feature is 
more obvious in Figs. [5] and IH] where we plot the differ- 
ential cross-section of pion and direct photon production 
at LHC as a function of rapidity at fixed transverse mo- 
menta Pt = 1,2,5 and 10 GeV within the GBW model. 
It is obvious that the invariant cross-sections have a peak 
at forward rapidity. However, compared to pions, the 
peak of the differential cross-section for direct photon 
production persists at larger pT- It seems that several 
mechanisms are at work here in different kinematic re- 
gions. Looking again at Fig. [7] it is obvious that in 
the case of direct photons when the saturation scale is 
smaller (the CGC model) at higher transverse momen- 
tum Pt = 2 GeV, the peak disappears and will be re- 
placed by a plateau. However, in the case of pion produc- 
tion, the peak is less pronounced even in the presence of a 
large saturation scale, see Figs. [7] and [S) Moreover, pho- 
tons are radiated by the electric current of the projectile 
quarks, which mostly stay in the fragmentation region of 



the beam, and tend to form a peak at forward rapidi- 
ties. However, at very large pT and 77, the kinematic 
limit pushes photon radiation to more central rapidities 
and the peak at forward rapidities will be replaced by a 
kind of plateau at central rapidities. At the same time, 
gluons are radiated via nonabelian mechanisms by the 
color current across the whole rapidity interval and tend 
to form a plateau at midrapidity. 

Another interesting difference between direct photon 
and hadron production is that direct photon production 
extends to higher rapidities for a fixed pT, see Figs. [8] 
and [HI This is more obvious in Fig. [TU] where we show 
the photon/pion ratio j/tt'^ as a function of pt at various 
rapidities within the GBW model and pp collisions. The 
ratio 7/7r" can be as big as 10 — 20 at very forward ra- 
pidities 77 ~ 8 — 7 at LHC energy. Note that suppression 
of hadrons at very forward rapidity also ensures signif- 
icant suppression of radiative decays of those hadrons. 
Therefore, direct photon production at forward rapidi- 
ties should be a rather clean signal. 

In Fig. Ill[ we show the ratio of photon/pion produc- 
tion as a function of rapidity in pp collision at LHC for 
various fixed pT within different saturation models. Di- 
rect photons can only be radiated from quarks, while 
hadrons can be produced by both gluons and quarks. At 
the LHC energy at midrapidity gluons dominate. There- 
fore the photon/pion ratio is significantly reduced to- 
ward midrapidity. However, at very forward rapidity, 
valence quarks become important and the photon/pion 
ratio rises. Moreover, at high pT again valence quarks 
becomes important and we have a sharp rise of the pho- 
ton/pion ratio, see Fig. [TT] A similar behavior has also 
been reported in a different approach [60| . 

In order to understand the relative importance of sat- 
uration effects at various rapidities, we employ the Semi- 
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FIG. 10: The ratio of photon/pion production in pp collisions 
at LHC as a function of pr calculated with the GBW color 
dipole model. 



10' 



10 



10 



^10-' 



10' 



10" 



; 1 1 1 1 1 1 1 1 1 


1 ' 1 ' 




: — GBW 
^ -- CGC 

■■■■ KLR-AdS/CFT 


14 TeV, pp 

p^=IOGeV 1 


P-,=5GeV : 
/ 1 




/ / / V 




/ y 

— .XV' 


/f y 

' f/ - 
// 






// 

GeV : 


. 'Tn^T^^' 






.1.1.1.1. 


I.I. 


. 1 . 



10 =— = 



1 2 3 4 5 6 7 



10-'^ - 



I ' I ' I ' I ' I ' I ' I ' I ' I 



Vr=14TeV,pp^Jt X X 



Semi-Sat, y =0.43 

*eff 

Semi-Sat, y =0.43+Diff 

'eff 




I ■ I ■ I ■ I ■ I ■ I ■ I ■ I ■ I 



2 4 6 8 10 12 14 16 18 20 
Pt [GeV] 

I ' I ' I ' I ' I ' I ' I ' I ' I ' 
Vr=14TeV, pp^yX J 

. Semi-Sat, Y ,,-=0.43 

'elt 

_ Semi-Sat, Y „=0.43-l-Diff 

'eff 




I . I . I . I 



2 4 6 8 10 12 14 16 18 20 
PT[GeV] 



FIG. 11: The ratio of photon/pion production in pp collisions 
at LHC as a function of rapidity r\ calculated with various 
color dipole models. 



FIG. 12: Pion (top) and direct photon (lower) spectra ob- 
tained from the Semi-Sat dipole model with two different ef- 
fective anomalous dimension 7e// at LHC and forward rapidi- 
ties in pp collisions. 



Sat model. In Fig. [T^ we show, the differential cross- 
section of pion and photon production at LHC, calcu- 
lated once with diffusion term and once without, i.e. 
7e// = 0.43. We recall that the Semi-Sat model in 
the presence of the diffusion term describes _F2-data at 
HERA, see Sec.VI-C. In Fig.[T2]we show that at forward 
rapidities, the diffusion term in the anomalous dimension 
is not important, since it gives similar results as with a 
fixed 7e// = 0.43. The preferred value of anomalous di- 
mension 1 — jeff — 0.57 at very forward rapidities is close 
to the one predicted from the BK equation (4^. This is 
more obvious in Fig. [T^] where we show the ratio of the 
two cross-sections for both pions and direct photons. It 



is well known that the saturation effects start being es- 
sential when the anomalous dimension reaches the value 
7cr = 1 — 7e// = 0.37 which is the case for forward ra- 
pidities (see Refs. 0,110, EH). This indicates that direct 
photon and hadron production at different rapidities at 
LHC are rather sensitive to saturation. 
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FIG. 13: The ratio of two invariant cross-sections obtained by 
taking two different effective anomalous dimension jef / in the 
Semi-Sat dipole model for pion and direct photon production 
at LHC energy (-^/s = 14 TeV) in pp collisions at various 
rapidities. 



VIII. CRONIN EFFECT AND NUCLEAR 
MODIFICATION FACTOR 



The nuclear modification (Cronin) factor Rpj^ is de- 
fined as ratio of pA to pp cross-sections normalized to 
the average number of binary nucleon collisions, 



R 



■pA 



dyd^PT 



V^b'inaryl dyd'-'pT 



(38) 



(Nbinary), the average number of geometrical binary col- 
lisions, is calculated according to the Glauber model [gl] 
for different centralities. 

Two very different mechanisms have been proposed to 
explain Cronin enhancement (or stippression) in pA colli- 
sions: a) initial-state effects [tI. [26l [63l - l65j due to a broad- 
ening of the parton transverse momentum in the initial- 
state. Here the fragmentation of hard partons is assumed 
to occur outside the cold medium, b) final-state effects 
Q due to the recombination of soft and shower partons 
in the final-state. 

In our approach, the Cronin effect originates from 
initial-state broadening of the transverse momentum of 
a projectile parton interacting coherently with a nuclear 
medium. The invariant cross-section of hadron and direct 
photon production in pA collisions can be obtained via 
the light-cone color-dipole factorization scheme defined 
in Eqs. (I^OIMl) . 



A. Gluon shadowing 

In the infinite momentum frame, the gluon clouds of 
nucleons which have the same impact parameter overlap 
at small Bjorken-x in the longitudinal direction. This 
allows gluons which originate from different nucleons to 
fuse, corresponding to a nonlinear term in the evolution 
equation which suppresses gluon production, and a preco- 
cious onset of the saturation effects for heavy nuclei. This 
is called gluon shadowing. The same effect, looks differ- 
ent in the rest frame of the nucleus, the gluon shadow- 
ing correction can be calculated as Landau-Pomeranchuk 
effect, namely the suppression of bremsstrahlung by in- 
terference of radiation from different scattering centers. 
This mechanism requires a sufficiently long coherence 
time for radiation, a condition equivalent to requiring 
a small Bjorken-x in the parton model. 

The question if glu on shadowing is an intrinsically 
leading twist effect [6a] or is suppressed by power of 
and is due to higher twist/high parton density effects 
[l|, is still debatable. There has been several shadow- 
ing models which consider only leading twist shadowing, 
e. g., including shadowing effects in the non-perturbative 
initial conditions which are then evolved with leading 
twist DGLAP equations 67). Modifications of this lead- 
ing twist picture to include Mueller- Qiu type non-linear 
contributions has been studied in Ref. [68|. 

In our approach, nuclear shadowing for gluons is cal- 
culated from shadowing of the \qqg) Fock component 
of a longitudinally polarized photon. Unlike transverse 
photons, all qq dipoles from longitudinal photons have 
size and the double-scattering term vanishes like 

l/Q*. The leading-twist contribution for the shadowing 
of the longitudinal photons arises, therefore, from the 
\qqg) Fock component. While the qq separation is of or- 
der the gluon can propagate relatively far from the 
qq-pair. After gluon radiation qq is in a color octet state, 
consequently the qqg system appears as gg dipole. The 
shadowing correction to the longitudinal cross-section is 
then directly related to gluon shadowing. The gluon 
shadowing ratio is defined as the ratio of the gluon den- 
sities in a nucleus and a nucleon nQ, ^ : 



RGix,Q^b) 



GA{x,Q\b) 
AGn{x,Q^) 



^^l%qg]{x,Q\h) 



(39) 

where lS.a'jf^[qqg] is the inelastic correction to the lon- 
gitudinal photoabsorption cross-section cr^^ due to the 
creation of a \qqg) Fock component. The details for the 
calculation of the suppression factor Rq can be found in 
Refs. [13, E^llSI- For a proton target, we have Rg = 1 
by construction. 

At high energy the qq dipole cross-section is also sub- 
ject to the multi-pomeron fusion effects in a nuclear 
medium. These effects are missed in the eikonal formulas 
Eqs. (j9ll5ll6ll7p where the variation of the transverse 
size of the qq Fock component while propagating and 
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interacting with a medium was not taken into account. 
Consequently higher Fock components were summed up 
without incorporating gluon shadowing. One should note 
that the multiple parton interactions that lead to gluon 
shadowing are also the source of gluon saturation. In 
order to avoid double counting, we calculate the nuclear 
shadowing effect within the same color-dipole formula- 
tion. The authors of Refs. [10, 69] have performed such 
a calculation by numerically solving the qq dipole evolu- 
tion equations in a medium by light-cone Green function 
techniques and confronted DIS data for nuclei. Follow- 
ing Refs. QjIiSllljIfiSllOl one can effectively incorporate 
gluon shadowing due to the nuclear medium by modify- 
ing the cross-section of the qq dipole interacting with a 
nucleus target at impact parameter b by the following 
replacement 



aqq{r, x) -> Raix, Q'^,b) x crgg(r, x), 



(40) 



in the exponent of Eqs. (j9ll5ll6ll7p . Therefore, by means 
of Rg and the dipole cross-section on a nucleon target, 
one can effectively define the qq dipole cross-section for 
a nucleus target by using Glauber theory, i.e. via simple 
eikonalization of the gg- nucleon cross-section modified by 
the suppression factor Rq. In this way, we relate the 
nuclear gluon shadowing to the gluon saturation which 
can be then read off from the constructed dipole-nucleus 
forward amplitude. However, the question if the parton 
saturation provides a precise microscopic understanding 
of shadowing is an open question and out of scope of this 
paper. One should also note that although the shadow- 
ing factor Rg improves the eikonal approximation, it is 
not apparently a solution of the non-linear BK evolution 
equation. 

In the CGC picture, the dipole-nucleus amplitude has 
the same functional form as the dipole-nucleon ampli- 
tude. The only difference is the saturation scale. The A- 
dependence of dipole-nucleus amplitude enters through 
the saturation scale Q\A^^^ where A is the effec- 

tive mass number of the nucleus in a given centrality and 
depends on the impact parameter. Our approach is dif- 
ferent but is not in contradiction with the CGC picture 
at the saturation boundary. In order to see this, let us as- 
sume that in spirit of the CGC picture one can write the 
forward dipole-nucleus amplitude in the following form 
(we use the GBW form for simplicity). 



(41) 



By comparing the above equation with Eq. (1151) and as- 
suming that there is no correlation between dipole am- 
plitude and the nuclear thickness, one can immediately 
read off the effective saturation scale in Eq. (|T5t close to 
the saturation boundary. 



Q^^(x,6) 



TA{h\ (42) 



2a^RG{x,Q')Qlj,{x)TA{h), (43) 



where in the second line we rely on the small-r ap- 
proximation of the dipole cross-section (valid for a large 
Ta(6)) and use the fact that Rg — t- 1 at r — )■ since 
q2 ^ Xjr^ — >. oo. Therefore, the square of saturation 
scale Q\j^ in our approach is approximately proportional 
to A^/"^ since Ta^) ^ A^l'^ ^ in agreement with the basic 
idea of saturation and the CGC picture fl|, |3, H, [zH • Let 
us repeat the above steps in a slightly different way. The 
dipole-nucleon cross-section at small dipole size r can be 
related to the gluon distribution xG(x^ Q^) in the nucleon 

[zl, 



^^,ir, x) = -a,{iy)xG{x, lyy. (44) 
By plugging the above expression into Eq. (|42l) we obtain, 

QlAix,b) = ^RG{x,Q')a,{l/r^)xG{x,l/r^)TAib), 

(45) 

where the typical value of dipole size can be related to the 
saturation scale Q^a ^ ^/r"^- This is remarkably similar 
to the saturation scale proposed by by Kharzeev, Levin 
and Nardi (KLN model) [tI], 

QIa{x, h) = ^-^a,{QlA)xG{x, Q^^)p^,,i(6), (46) 

where for pA collisions the density of participants is 
Pparti^) — TA{b)- Note that Eq. p6)) gives the saturation 
scale for gluons and it is different from the saturation 
scale for quarks Eq. (|45|) by a Casimir factor 9/4. The 
KLN model Eq. pS)) gives a good description of hadron 
multiplicities in heavy ion collisions at RHIC (t^]. The 
main difference between our model Eq. (j45|) and the KLN 
model Eq. is the shadowing factor Rg{x, Q^) which 
takes into account approximately multi-pomeron fusion 
effects in a nuclear medium beyond the eikonal approxi- 
mation, see also Ref. |74|. 

In the limit of strong shadowing at very small Q^, 
the gluon ratio Eq. ([M]) has a simple form Rg ~ 
ttRa I i^'^ef f ) where a^f j is the effective cross-section 
responsible for shadowing and Ra is the nuclear radius. 
Therefore, in our approach, deep inside saturation region, 
we have Rg ^/TA(b) and consequently the saturation 
scale QsA becomes independent of A. This behavior has 
been also predicted based on more sophisticated mod- 
els indicating that the parton wave functions of different 
nuclei become universal at high energies limit [75j, |76| . 
Nevertheless, we expect that our approach based on an 
improved eikonal approximation will not be reliable at 
such an extreme limit and we use our formulation only 
at midrapidity for pA collisions at RHIC and LHC ener- 
gies. 

As we argued above, in principle one may construct 
the dipole-nucleus amplitude via the dipole-nucleon am- 
plitude supplemented with the A-dependent saturation 
scale. However, it is not a priori obvious whether such a 
model with parameters fitted to the available DIS data 
on proton target is also able to describe the DIS data 
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FIG. 14: Nuclear modification factor RpA for pion produc- 
tion at RHIC energy in minimum bias proton-gold collisions 
at midrapidity within the long-coherence length scheme pre- 
sented in this paper. Note that at RHIC energy we are in 
the transition region between the short and long-coherence 
length limit. For comparison, we also show the Cronin curve 
obtained in the short-coherence length scheme [s^. The 
Cronin ratio for gluon production in the long-coherence length 
scheme is also shown. The GBW saturation model is used for 
all curves. The experimental data are from [53l ]. 



on nucleus target at small-x without having to change 
the parameters of the model (the issue of sensitivity of 
the model parameters obtained from a fit to data in 
analysis), see Ref. [77[. In our approach we use the same 
dipole-nucleus cross-section which gives a good descrip- 
tion of HERA data to calculate the cross-section in pA 
reactions. We stress again that the shadowing factor 
Rg{x,Q^) is not a free parameter in our formalism but 
it is calculated via Eq. ([39t . Such a shadowing factor 
is needed in order to describe the DIS data off nuclei 
0, HO, m, lei, [zO]. Therefore the suppression obtained 
as a result of the inclusion of the shadowing factor Rq 
(which depends on kinematics) is not arbitrary. 



B. Numerical results for pA collisions 

For the calculation of cross-sections for pA collisions, 
we use the same PDFs and FFs as for pp collisions. 
Furthermore, we use a Woods-Saxon nuclear profile for 
TA{b). We again stress that similar to the calculation for 
pp collisions, here again we have no free parameters to 
adjust. In Fig. [HI we show RdAu for 7r° production at 
RHIC in minimum bias proton-gold collisions. The ex- 
perimental data in Fig. [HI are from PHENIX [53]. As we 
already mentioned, for RHIC energy at midrapidity and 
moderate p-r, the coherence length defined via Eq. ([2]) is 
about Ic ~ 5— 6 fm which is comparable to the nuclear ra- 
dius. Therefore, we are in the transition region between 



FIG. 15: Nuclear modification factor RpA for pion production 
at the LHC energy ^/s — 5.5 TeV at midrapidity in minimum 
bias proton-lead collisions within the CGC and GBW color 
dipole models. We show results with and without inclusion 
of the nuclear gluon shadowing factor Rg{x, Q^,b). 



the regimes of long and short coherence length. Calcula- 
tions in such a region are most complicated. In Fig. I14[ 
we show the theoretical curves calculated in the two ex- 
treme cases of short- and long-coherence length. The 
curve for the short-coherence length in Fig. [T3] is based 
on an improved pQCD calculation taken from Ref. [ssf . 
We used the AKK08 for FFs, MSTW2008 for PDFs and 
the GBW model for the color dipole cross-section. One 
should also note that the color dipole cross-section is fit- 
ted to the DIS data for X2 < 0.01. Therefore, our results 
at high pt for RHIC are less reliable. At RHIC and LHC 
energies at midrapidity, gluons are mostly responsible for 
pion production. We also show in Fig. ll4l the Cronin ratio 
for gluon production in proton-gold collisions. It is seen 
that fragmentation processes distort the gluonic Cronin 
enhancement and shift the Cronin peak to a lower px- 

In Fig. [13 we show our prediction for the nuclear mod- 
ification factor RpA for tt" production at LHC at midra- 
pidity in minimum bias pA collisions within two very 
different saturation models, namely GBW and CGC. We 
also show the effect of nuclear gluon shadowing. It is seen 
that the Cronin enhancement will be replaced with mod- 
erate suppression in all saturation models considered in 
this paper due to nuclear gluon shadowing. It is obvious 
that a qq-proton dipole model with a bigger saturation 
scale leads to a larger Cronin enhancement and works 
against the nuclear shadowing suppression. This effect 
has also been shown in Ref. J7^]. Note that the source 
of both saturation and shadowing is parton multiple in- 
teraction. However, a larger saturation scale leads to a 
stronger broadening of transverse momentum of the pro- 
jectile partons and consequently it works against shad- 
owing. This is more obvious in Fig. (jl6p (uppor panel) 
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FIG. 16: Same as Fig. [15] for gluon production. In all curves 
in the upper panel the nuclear gluon shadowing factor Re 
(defined via Eq. (I39p ) is incorporated. Lower panel: gluon 
shadowing effects at LHC for the GBW model. 



FIG. 17: Same as Fig. [15] for direct photon production. For 
comparison, we also show the results with and without in- 
clusion of nuclear gluon shadowing effects introduced via 
Eq. ((391). 



where we plotted the Cronin ratio for gluons produc- 
tion at the LHC energy within various saturation cofor 
dipole models. In Fig. (lower panel) we show effect 
of nuclear gluon shadowing within the GBW color dipole 
model. It is seen that both shadowing and saturation 
effects are important at LHC in pA collisions and give 
rise to a rather sizable effect in the nuclear modification 
factor RpA- 

In Fig. [iTl We show our prediction for the nuclear 
modification factor i?^^ for direct photon production at 
LHC at midrapidity in minimum bias pA collisions for 
two models with different saturation scale. In order to 
demonstrate the importance of nuclear gluon shadow- 
ing effects, we have also plotted the curves without nu- 
clear gluon shadowing. In comparison to pion produc- 
tion, the Cronin enhancement for direct photon produc- 
tion seems stronger and survives within the GBW color- 
dipole model which has a bigger saturation scale, even af- 
ter the inclusion of nuclear gluon shadowing suppression 
effects. Similar to pion production, the Cronin enhance- 
ment for direct photon production is bigger in a model 
with a larger saturation scale. Within the CGC model 
both pion and direct photon enhancement at RHIC will 
be replaced by suppression at LHC. 

In a similar approach, Kopeliovich et al. have 
shown that the Cronin enhancement will survive at LHC 
at midrapidity though reduced compared to RHIC. Here, 
our finding is different. This is due to the fact that we im- 
proved the earlier calculation in several ways including: 
using updated PDFs, FFs and color-dipole cross-sections, 
incorporating gluon radiation from the quark projectile 
(i.e. qN — )• qgX, see Eq. (I^ 'l. us ing recently upgraded 
shadowing suppression factor Rg [TOj and using a run- 
ning strong coupling. Nevertheless, in both approaches 



the Cronin ratio RpA is still very small, less than 20% at 
midrapidity at LHC. Kharzeev et al. [64] have shown a 
marked suppression for pions at midrapidity at LHC in 
pA collisions based on the CGC scenario. This suppres- 
sion is stronger than our prediction. Certainly, LHC data 
should be able to decide between the different approaches 
and scenarios. 

Notice that our prescription for both hadron and pho- 
ton production in pA collisions is less reliable at very 
large pr and also forward rapidities. This is due to the 
fact that at large xp (i.e. Xi — )■ 1) one should properly 
incorporate energy conservation since it puts an impor- 
tant constraint on particle production Q . Nevertheless, 
we expect this effect to be negligible in our kinematical 
region of interest. Note also that the energy loss effects 
are subject to xi-scaling and less important for high- 
energy pA collisions at moderate pr, although it might 
be important at lower energies 0, [l^ . A more detailed 
study of the Cronin effect for direct-photon production 
at RHIC and LHC at forward rapidities and high pT will 
be presented elsewhere [7§|. 



IX. CONCLUSIONS AND OUTLOOK 

In this paper, we investigated pion and direct pho- 
ton production within a unified color-dipole approach at 
high-energy pp and pA collisions and provided various 
predictions for the upcoming LHC experiments. The re- 
sults of this paper can be summarized as follows: 

• Both hadron and direct photon production strongly 
depend on the value of the anomalous dimension 
7e// and are sensitive to gluon saturation effects 
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at forward rapidities at LHC (^s = 14 TeV). The 
difference between various saturation model predic- 
tions can be about a factor 2-^3. Note that all sat- 
uration models employed here are fitted to HERA 
data. 

• We showed that the ratio of photon/pion produc- 
tion at LHC {^/s — 14 TeV) at very forward ra- 
pidities in pp collisions can be as big as 10 — 20. 
Therefore, direct photons at very forward rapidi- 
ties should be a rather clean observable and provide 
a sensitive probe for saturation effects and small-x 
physics in general. 

• We showed that the rapidity distribution of pi- 
ous and direct photons exhibit some peculiar en- 
hancement at forward rapidities which is more pro- 
nounced in the case of photon production. This 
peak is enhanced in models with a larger satura- 
tion scale at lower px ■ 

• We investigated the relationship between satura- 
tion and shadowing effects in pA collisions at LHC 
for both direct photon and hadron production. We 
studied the role of initial-state broadening of the 
transverse momentum distribution of a projectile 
parton propagating and interacting coherently with 
a nuclear medium. We showed that a larger satu- 
ration scale leads to a stronger transverse momen- 



tum broadening of the projectile partons and conse- 
quently works against the nuclear gluon shadowing 
suppression effects. Our results show that the nu- 
clear modification factor RpA at LHC is sensitive to 
both saturation and nuclear shadowing effects and 
it seems that a subtle cancellation between these 
two effects leads to a rather small Cronin ratio RpA- 
We showed that the tt'^ and direct photon 7 Cronin 
ratio RpA at the LHC is less than 1 within the CGC 
color dipole model. However, in the case of direct 
photon production in pA collisions, the Cronin en- 
hancement can survive at the LHC energy within 
the GBW color-dipole model which has a larger 
saturation scale. 
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